Analytical Study of Mode Coupling in Hybrid Inflation 
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We provide an analytical study of the coupling of short and long wavelength fluctuation modes 
during the initial phase of reheating in two field models like hybrid inflation. In these models, there 
is - at linear order in perturbation theory - an instability in the entropy modes of cosmological 
perturbations which, if not cut off, could lead to curvature fluctuations which exceed the current 
observational values. Here, we demonstrate that the back-reaction of short wavelength fluctuations 
is too weak to lead to a truncation of the instability for the long wavelength modes on time scales 
comparable to the typical instability time scale of the long wavelength entropy modes. Hence, unless 
there are other mechanisms which truncate the instability, then in models such as hybrid inflation 
the curvature perturbations produced during reheating on scales of current observational interest 
may be very important. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The inflationary scenario [l[ is the current paradigm 
of early universe cosmology. It addresses several con- 
ceptual problems which Standard Big Bang cosmology 
was unable to solve, and - perhaps more importantly - 
provides a causal mechanism to generate the primordial 
cosmological fluctuations which lead to the large-scale 
structure and cosmic microwave background anisotropies 
which are currently observed @ • Reheating after the pe- 
riod of inflation is a key aspect of inflationary cosmology 
(see the recent review of [3]). During reheating the en- 
ergy is transformed from the inflaton, the scalar matter 
field which is responsible for providing the inflationary 
expansion of space, to the matter we see today. Without 
reheating, the inflationary universe would leave behind a 
large universe empty of any matter, in obvious contradic- 
tion with observations. Reheating requires some coupling 
between the inflaton and regular matter. 

The initial studies of reheating were based on first or- 
der perturbation theory applied to the inflaton decay 0] ■ 
This approach, however, fails to take into account the 
coherent nature of the inflaton field. At the end of the 
period of inflation, the inflaton begins coherent homoge- 
neous oscillations about the ground state of its potential. 
In Q (see also [1, 0] ) it was realized that this gives rise 
to a potential parametric instability in any matter field 
which couples to the inflaton. This instability can take 
various forms depending on the nature of the coupling 
between the inflaton and the matter field. Rather gener- 
ically, the process leads to a rapid energy transfer from 
the inflaton to the matter fields [4l[. In a wide class of 
models, the resonance is of "broad resonance" type 
and affects all fluctuations with wavenumbers smaller 
than a characteristic mass set by the particle physics La- 
grangian. In some cases - in particular in the hybrid 
inflation model which we study in this paper - there is 
a tachyonic instability in the matter sector at the end 
of inflation, in which case the reheating instability is of 
"tachyonic resonance" Q type and hence extremely effi- 
cient. 



The inflaton field obviously couples to gravity, and 
therefore, as first conjectured in [10], it is possible that 
its oscillations will lead to a parametric resonance insta- 
bility of the cosmological fluctuation modes. Due to the 
exponential growth of the causal horizon during the infla- 
tionary phase with constant Hubble radius, even super- 
Hubble modes can be causally affected by this instability 
11]. This will inherently cause a long wavelength curva- 
ture mode to develop which potentially exhibits a tachy- 
onic resonance. These large scale modes are of prime 
importance, since they correspond to scales observable 
today. If an inflationary model fails to control this possi- 
ble instability of the curvature mode, the observed small- 
ness of the magnitude of curvature perturbations today 
might rule out this model. 

If matter consists of a single scalar field, then - as stud- 
ied in [TlT | and 12 1 - there is no parametric instability of 
super-Hubble scale modes. However, if entropy fluctua- 
tions are present, then an instability of entropy modes 
can occur which will lead to rapid growth of the curva- 
ture fluctuations. For this instability to be effective, the 
entropy mode cannot have been exponentially redshiftcd 
during the period of inflation. Models which satisfy the 
conditions for parametric instability of the metric fluctu- 
ations were first discussed in [l3[ and [l4[ . 

If the duration of the parametric instability of long 
wavelength curvature fluctuations were long, then the 
curvature fluctuations induced by the resonant entropy 
modes could become larger than the primordial curvature 
perturbations and - in fact - larger than the observed am- 
plitude of the inhomogeneities. Thus, one might poten- 
tially be able to use the parametric resonant instability 
of the entropy modes of the cosmological fluctuations to 
rule out large classes of multi-field inflationary models. 

One class of inflationary models which at linear or- 
der in perturbation theory leads to an instability of the 
entropy fluctuations is hybrid inflation [15j. In this sce- 
nario, two scalar fields <f> and tp, are involved in the infla- 
tionary process, the slowly rolling inflaton field <p an( l the 
"waterfall" field ip. The key point of this class of mod- 
els is that while <j> is slowly rolling during inflation, the 
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energy density of the Universe is dominated by the po- 
tential of ip. Over the past few years, a lot of interest has 
been devoted to these models, mainly since they provide 
a framework for the realization of inflation in the context 
of both supersymmetry [l6[ and string theory (see [l?} 
for reviews and [l8| for an original reference). Among 
the promising approaches stand the D-brane / antibrane 
inflation models (e.g. the D3/D7 brane inflation model 
19] and the KKLMMT model [IS]). Though these mod- 
els may resolve some of the conceptual problems from 
which simple single scalar-field driven inflation models 
suffer, the large number of light moduli fields present 
in string compactifications can give rise to entropy fluc- 
tuation modes, which could in turn enter a parametric 
resonance phase at early stages of reheating [l3|, 13 • 

Hybrid inflation models are characterised by symmetry 
breaking along the ip direction: The symmetry ip — > —-ijj 
is spontaneously broken for field values <p smaller than a 
critical value <p c . Inflation takes place while <fi is slowly 
rolling at field values <p > 4>c- Once the inflaton crosses 
the critical value, the effective potential for ip develops 
a tachyonic instability which triggers the rapid rolling of 
ip towards one of the ground states of its potential. This 
causes inflation to end. 

Reheating in hybrid inflation models on a fixed Fried- 
mann background cosmology (i.e. no cosmological per- 
turbations) was studied in detail in [2l|,[22[ (see also [Hj]). 
By means of numerical simulations it was observed that 
within a very short time, non-linearities on a length scale 
given by the mass of the waterfall field develop and dom- 
inate the subsequent stages of the reheating process. It 
is important to study how large-scale metric perturba- 
tions evolve during reheating in these models. In pre- 
vious work, linear evolution of the entropy and curva- 
ture modes was studied 0, HH [HJ , and it was shown 
that it is possible that an important entropy perturba- 
tion mode develops. In these works it was assumed that 
non-linear effects, e.g. the back-reaction of small-scale 
fluctuations, would unlikely have a dominant effect on 
large-scale fluctuation. Due to the large range of scales 
(the small wavelength modes we are interested in have 
wavelengths comparable to the inverse Hubble radius at 
the end of inflation whereas the wavelengths of modes 
we are interested in for current cosmological observations 
are of the order of 1mm, assuming that the scale of in- 
flation is about 10 15 GeV) numerical studies do not have 
the dynamical range to study this question. Instead, an 
analytical understanding of the dynamics is required. 

Our work demonstrates, by means of an analytical 
analysis, that the back-reaction of short wavelength 
modes with random phases is too weak to truncate the 
instability of the long wavelength entropy modes. We 
have studied the effects of short wavelength fluctuations 
both on the background inflaton and on long wavelength 
modes of the inflaton and waterfall fields. Most impor- 
tantly, we have shown that the back-reaction on large 
scale fluctuations of the waterfall field is too weak to shut 
off the resonance of this entropy mode one the relevant 



time scale of the instability. This implies that hybrid 
inflation models of the type analyzed in this paper may 
suffer from a potential "entropy mode problem" , unless 
there are other ways of truncating the resonance (e.g. 
the exponential decrease in the initial value of the fluctu- 
ations modes during the period of inflation - see e.g. [25j 
for a discussion) [43J . 

Our study of the back- reaction of fluctuations in hybrid 
inflation models is not the first analytical study. For an 
analytical study of the back- reaction effect of fluctuations 
on the background u sing very different techniques than 
the ones we use see [31(. For a study of the effects of 
metric fluctuations on the observable local expansion rate 
of the universe due to fluctuations see [l^l ■ 

The structure of this paper is as follows. In Section 2, 
we review the class of hybrid inflation models which are 
studied here, and introduce the perturbative approach 
employed to study non-linear interactions that arise dur- 
ing preheating. In Section 3, we solve the background 
theory, and in Section 4 we review the evolution of quan- 
tum fluctuations to linear order. In Section 5, we analyze 
the generation of long wavelength second order pertur- 
bation modes sourced by shorter wavelength first order 
fluctuations (summing over the contribution from all fre- 
quencies and assuming random phases), for both the in- 
flaton and the waterfall field. This is the leading order 
back-reaction effect between short and long wavelength 
modes. We also study the back-reaction of the fluctua- 
tions on the background. We find that small-scale first- 
order modes with random phases have a contribution on 
the evolution of large-scale modes which does not start 
to dominate until long after the instability of the first 
order long wavelength modes has had a chance to de- 
velop. Hence, in spite of their large phase space, the 
non-linear effects of the short wavelength perturbations 
do not become dominant for the evolution of long wave- 
length modes during the early stages of preheating. 



II. THE MODEL AND PERTURBATIVE 
APPROACH 

We focus on the category of hybrid inflation models 
with Lagrangian density for matter given by: 

C m {^) = Id^d^+^d'rjj (1) 

The equations of motion form a coupled system of partial 
differential equations 

4) + 3H0- \\7 2 (j) = - (m 2 + g 2 ip 2 ) <f> (2) 
i> + 3Hip- \v 2 ip - - (A(V> 2 -v 2 )+g 2 (f> 2 )ip.(3) 
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The Hubble parameter H (t) at time t is given by: 



, a\ 8ttG 
H{t) = [a) = ~ >' 



where p is the energy density, and is given by: 



(4) 



(5) 

where the potential V is given in the second line of ([1}. 
The turnover value of 4> at which ip develops a tachyonic 
instability will then be 



vX 1 / 2 
9 



(G) 



To study this model we work in discrete Fourier space 
(discrete because of a finite volume cutoff) and expand 
to second order about a homogeneous and isotropic cos- 
mological background. The expansion parameter e is the 
amplitude of the linear fluctuations. Our goal is to study 
how first order perturbations of high k modes (with wave- 
lengths comparable to the Hubble length at the end of 
inflation or to the wavelength associated with the mass of 
the waterfall field, whichever is smaller) influence the low 
k modes (modes which affect cosmological observations 
today which correspond to a scale of roughly ~ 1 mm 
at the end of inflation) at second order. We want to es- 
timate the time interval At it will take before this back- 
reaction effect becomes dominant, and we also want to 
see whether back-reaction effects on modes with wave- 
length of the order Iq decreases with the wavelength. 

The expansion of the fields to second order in e is: 

<t>(x,t) = (t) {0} {t) + e5(t) {1) (x,t) + e 2 S(t) {2) (x,t) (7) 
ip(x,t) = + eS^ 1 \x,t)+e 2 Sij {2} (x,t). (8) 

Since we will eventually be evaluating mode sums nu- 
merically, it is useful to work with real Fourier modes. 
Hence, we can expand the first and second order field 
perturbations as follows: 



Here, 2L is the size of the finite one-dimensional box 
inside of which we are performing the discrete Fourier 
expansion. Physically, we need to take L larger (but not 
much larger) than the largest scale of the problem we 
study; hence we fix it to be ~ 1mm = 6 x I0 31 l p . This ef- 
fectively imposes a cutoff on the largest scale studied. For 
clarity, we only write explicitly the one-dimensional ex- 
pansions, but generalisation to three-dimensional Fourier 
series is straightforward and will not modify in a crucial 
way the form of the obtained solutions, unless otherwise 
mentioned. 



III. ZEROTH ORDER EXPANSION 

Inserting the above ansatz into the equations of motion 
of the system and expanding to zcroth order in e, © and 
d2]) reduce to 

(o) (i)+3#0 (o) (i) = -m 2 (f> (0 \t) (9) 
V> (0) (i) = 0. (10) 

For values of \<f>\ smaller than </> c , there is an instability 
of the background solution for ip. Because of this insta- 
bility, the ip field grows fast on the scale of a Hubble 
expansion time. Hence, we expect the Hubble damping 
term to be negligible, and thus we can set H sa 0. This 
amounts to setting a(t) to a constant (which we pick to be 
1) and a(t) = 0. The linear equation for <f> then becomes 
that of a harmonic oscillator and has the solution 



(o) (i) = A (0) cos(mt) + B (0) sin(mt) , 



11) 



where and B^ are constants depending on the ini- 
tial conditions on <^°) (t = 0) and its derivative. Since in 
our case, we are interrested in the end of the inflationary 

r and we also want d> to be 



era, we want <f>(°'(t = 0) 
initially slowly rolling, i.e. 
set B^ to zero and obtain: 
0(°)(t) = 



0<°>(t = Q) 

4> c cos(mt) 
0. 



< 1. Thus, we 

(12) 
(13) 



COO 
EK 1)8 (*)-(f) +?i " )c(<)cos (f) 
n=0 



<M 2) = 



Sip 



(2) 



/ oo 

E 

\n=0 



/ OO 

\n=0 



f oo 
\n=0 



4? c (t)cos (- 



+ ^(t)cos( 



V L 



+ ^ c (t)cos( 



V L 



IV. FIRST ORDER EXPANSION 
A. Equations 

Now, going back to the system @ and © and insert- 
ing the ansatz (J7J and (|5J|, we expand and keep terms of 
first order in e: 

54> {1 \x,t) + ZH5ft l \x,t) - ^V 2 ^ (1) 0M) 

= -m 2 5<j) {1) {x,t) (14) 
5ii>W(x,t) + 3HSip ( - 1) (x,t) ~ -^V 2 <M> (1) (M) (15) 

= Xv 2 5^ 1 \x,t) - g 2 ^°\x,t)) 2 (S^ 1 ^) . 
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Inserting the explicit form of the first order pertur- 
bations, we make use of the orthogonality relations for 
trigonometric functions to convert (fl"4"|) and (fT5)l to dis- 
crete Fourier space. 

Doing so, we obtain the following set of differential 
equations for the first order correction to each Fourier 
mode: 



The equation can be reduced to the Mathieu equation 
by performing the transformation z = mt and using the 
identity cos 2 (z) = 1/2 (1 + cos(2z)): 



€ (1)s ' c + 



3H<p. 



= 



(l)s,c 



(l)s,c 



/nir\- 
\oLJ 



/nir\ • 
\aL) 



A, 2 +. 9 2 (>>) UO 



= 0. 



Replacing the wavenumber n by its vectorial expression ii 
yields, without any other modification, the generalisation 
of the 1+1-dimensional equations to the 3+1-dimensional 
setting. 



B. 



Solutions 



Again setting H w 0, the first order <f> equations be- 
come that of harmonic oscillators, and so can be solved 
easily: 



(l)s,c 



^COS 



© +mH . 



+£i 1)s 



) + m 2 t\ . (li 



Let us now have a closer look at the characteristic pa- 
rameter values. We have in mind a hybrid inflation model 
stemming from string scale physics. Hence, the value of 
v will be taken to be v = 10~ 2 in Planck units. We 
choose coupling constants g = 10 -2 and A = 10~ 4 , so 
that 4> c = 10 -2 (again in Planck units). We will consider 
the range m e [10 -7 , 10 -5 ] (in Planck units). For a value 
of m at the upper end of this interval, the hybrid infla- 
tion model will result in cosmological fluctuations of the 
observed order of magnitude (see e.g. [34j for a review of 
the theory of cosmological perturbations and [35| for an 
introductory overview) . We explore a range of masses m 
in order to study how the strength of our back-reaction 
effect depends on the model parameters. 

Moreover, we are obviously interested in modes whose 
wavelength is larger than m _1 and larger than the Hub- 
ble radius. For such modes (^) <C m. Therefore, (fT8|) 
can be approximated by: 



A^ )s ' c cos (mt) + B^ s ' c sin (mt) , (19) 

which describes stable harmonic oscillation. No instabil- 
ity is manifest in the <fi held. 

On the other hand, the first order ip equations become: 



&""(<)- 



(Hi) 



g 2 4> 2 c cos 2 (mt) 



^ )s ' c (t) = o. 
(20) 



/ 7J7T \ ' 

\amLJ 



m 2 



2a2 



9 4> \ 
2m 2 



(1 + cos(2z)) 



(1)8,0 



(21) 

where the prime refers to a derivation with respect to z. 
Defining: 



(16) 



',(1)8,0 



(17) 



4m 2 ' 

/ 717T \ ' 

\amL ) 



X<P_ 

m 2 



2q 



(22) 
(23) 



we indeed recover the canonical form of the Mathieu 
equation (see e.g. (33|) 

<« s ' c (i) + K - 2qcos(2z)] ^ s ' c (t) = . (24) 

The paramter q is called the "Floquet exponent" , and we 
will call oj n the "square frequency" . 

For the parameter values we are using, the value of q 
is much larger than 1. Hence, we are in the parameter 
region of either "tachyonic resonance" (if the tachyonic 
term in the expression for uj n dominates over the third 
term, or "broad parametric resonance" if the third term 
dominates (the first term in ui n is negligible for the modes 
we are interested in). In either case, all of the infrared 
modes which we study here will experience an exponen- 
tial instability with a growth rate characterized by the 
Floquet exponent. To first order in perturbation theory, 
every ip mode evolves in an independent way and there 
is no interaction between different modes. In particular, 
a mode that is not initially excited will not grow to first 
order at any later time (obviously, we expect quantum 
vacuum fluctuations on all scales to seed the instability) . 
Inserting the expression © for (j) c , we find that for all 
values of n of interest to us 



/ mr \ 2 
\amLJ 



A^ 2 
m 2 

9^1 



9 2 € 



2m 2 



2m 2 



(independently of the value of m); which means that uj n 
is always negative. Hence, we conclude that all modes 
we are interested in undergo tachyonic parametric reso- 
nance. 

The solution to the second order differential equation 
(|24p can be written in terms of two linearly indepen- 
dent solutions, the so-called Mathieu functions M at hC 
and M ath S: 

^ )s - c (z) = Ci 1 ^ c M ath C(oj n ,q,z)+Dl 1 ^ c M ath S(u n ,q,z), 

(25) 

where the C's and D's are the coefficients. 

Note that an important property of the solution to 
(|24l) for any choice of the parameters is the existence of 
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an exponential instability of parametric resonance type 
for certain ranges of the parameter value uj n . For values 

of uj n falling within the instability range, (j>n c (t) increa- 



we obtain 



220- 
200- 




FIG. 1: Analytical solution to (f24|) in 

red, cosh (\ _ cos(ms))) in green, and 
exp [^^(1 — cos(ms))~\ in yellow, all on a logarithmic 
scale, as a function of z = ms; for m = 10~ 6 and for the 
reference parameters cited in the text. The cosh function is 
a better approximation to the MathieuC function, and it is 
an upper bound for z < n. 

ses exponentially. In the cases of tachyonic and broad 
parametric resonance 



4! )s ' c (t) cx exp(nz) 



(26) 



for some constant fi depending on q. In the case of broad 
parametric resonance [8j , the parameter /x is of the order 
(but slightly smaller) than 1. 

We can find an approximate solution of (|24p valid over 
the first half of the period of under the assumption 
that the initial value is 1 (this is a normalization) and 
that the initial velocity of the mode vanishes (which, as 
discussed in the next subsection, is a good assumption 
for the modes we are interested in). Our approximate 
solution in fact gives an upper bound on the value of the 
mode function (which is within a factor two of the exact 
solution) and is given by 

M a thC\--— 7 ,— — 5-, z) < cosh I (1- costs)) 

\ 2m z Am* / \ m 

. (27) 

This result is found by approximating the equation of 
motion as 



) sin zi\)-\ sinzsinh — cosz 



But ~ 100 from our choice of parameters, and 

sinha; < coshx. So we have ib > ib at least from z = 
to z — 7r (in which region cosz > 0). 

C. Initial Conditions 

Initial conditions on the first order fluctuation modes 
are given by the quantum fluctuations at the end of in- 
flation. These modes begin on sub-Hubble scales at the 
beginning of the inflationary phase in their quantum vac- 
uum state. As reviewed e.g. in [34], [HI, the fluctuations 
freeze out and undergo squeezing once the wavelength ex- 
its the Hubble radius. The squeezing implies that the ve- 
locity of the mode functions will redshift. For purely no- 
tational simplicity we have chosen to excite only the (j>^ c 
and ij)V-> c modes and to set the ^ < - 1 - )s , tp^ s modes to zero. 
This implies that we are taking correlated phases for the 
first order modes. At the end of the calculation we will 
restore the randomness of the phases and comment on the 
effect that this has on the strength of the back-reaction. 
Because of the squeezing of the super-Hubble modes dis- 
cussed above, we start <jP-'° and ip^ c with zero velocity, 
that is to say, we set Bn° — Dn = V rt, while A^ c 
and Cn are determined by the Bunch-Davies state. 

When numerically computing mode sums later on in 
this article, it is important to know the initial amplitude 
of the mode functions in discrete momentum space. In 
continuous Fourier space, these amplitudes are given by 
(in d spatial dimensions): 

(1)c (k) = fc~ 1/2 , (30) 
where we are using the Fourier decomposition in the form 

d d k 



(i) 



(s,0) = me 



(2ir) d 



i7rlk:22(l) t 



{Tk)V 



1/2 



. (31) 



where V is the spatial volume. We want to match this 
set of initial conditions with our discrete Fourier series: 



'(=,0) = me 



E 

.11=0 



^4 1)C («) 



(32) 



(1 — cos 2 z)tp 



sin 2 (z)V>, (28) 



and imposing the initial conditions mentioned above. In- 
deed, if we set 



ib = cosh (1 

1 m 



cos z 



(29) 



We know TkL = 7rin and AkL = n, and thus want to 
relate ^' c (m) to ^^^(k) in terms of m. To do so, we 
make use of the identity (recalling that V = (2L) d ) 



me 



1 oo 

-T 



V 



me 



E 







' f d d k ' 




-> me 




[J (2n) d \ 



(33) 
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in the continuum limit in d dimensions. In this limit, the 
discrete expansion (|52")l needs to converge to (|3"Tj) . i.e.: 



Die 

me 
me 



Ak\ d 
~2^J 



E ( ^ ) ^ )e w 



(i#^ 1)c (^ 



2tt 

rl/2 



27T 



(34) 



where the last step expresses our requirement of conver- 
gence. Hence, for the initial values of the discrete Fourier 
modes we find the relation: 



^ 1)c (fc,o) 



^) $>>(») (21) 



d/2 



1 

2L 



1 



1 



(2tt) 1 /2 



(35) 



However, even though initial velocities of all first order 
modes must be zero due to freezing outside the Hubble 
radius, their relative phases must, in general, be random. 
Thus, the initial conditions must include these random 
phases. When measuring a first order mode, the expec- 
tation value of its amplitude in absolute value needs to 
be taken, which will divide the obtained amplitude by 2, 
while for further calculations, keeping the phase general 
as a random variable will be required. 

Bringing everything together, we thus write the solu- 
tion for the first order <fi and tjj modes: 



cos(6>„) 




J 



V. SECOND ORDER EXPANSION 

A. Equations 

Going back to the system © and ([3]) and again in- 
serting the ansatz ([7]) and (|8]), we expand and now keep 
terms of second order in e. We obtain 

54> {2) [x,t) + 3H8^ 2) (x, t) - \ V 2 (50 (2) (x, t) (38) 



^m 2 S^ 2 \x,t) - g 2 (81>W(jc,t) 



and 



8^ 2) (a, t) + 3H6iP (2) (:r, t) - \ V 2 <5^ (2) (x, t) (39) 



A. 2 - 5 2 (0 (o) ) : 



S^ 2 \x,t) 



2g 2 5^(x,t)5^(x,t)4>^ 



Inserting the explicit form of the first and second order 
perturbations, we make use of the orthogonality relations 
for trigonometric functions to convert (|3"8"f and (|3T)f to 
discrete Fourier space. However, this time the process is 
slightly non-trivial due to the presence of the interaction 
terms at this order in perturbation theory which give rise 
to mode mixing. Indeed, the last terms in equations 
and (|39)l describe how the growth of first order perturba- 
tions will source second order pertubations. They involve 
products of modes, which requires the use of trigonomet- 
ric identities to split these terms in a way that allows the 
use of the canonical orthogonality conditions for sines 
and cosines. After some algebra, this yelds the follow- 
ing set of differential equations for the second order cor- 
rection to each Fourier mode. In one spatial dimension 
(recalling that (f>Q S (t) — ^^(t) = by definition, and 

0o^ C (*) = V'o^W = because these modes are part of 
the background) we obtain the following results: 



For n > 1, the equations describing the back- reaction of the fluctuation modes on the perturbations themselves take 
the form 

2 00 

. $»>•(*) + 3H^(t) + 4? S (t) = ~m 2 ^(t) g 2 (*<°>(t)) £ [^(t) (^(t) _ ^ (t) )] (40 ) 

3=1 

. ^)«(t) + 3F^)«(t) + 0%! 2)C W = -m^Ht) 

2 00 

- y (^ (0) w) E W )s ® 0l-«(*) +4+n(t)) +4 )c ® ( 41 ) 



3 = 1 



7 



OO 

9 2 (* (0) (*)) E [(C»|(*> - ^nW) >- (*) + W£°n|(*) - >- (*)] (42) 



fc=l 



OO 

3 2 (> (0) (*)) E K 1)s w +^ 1)c w (^w +^n(*))] . (43) 



fc=i 



while for n = 0, that is, for the back-reaction on the background fields, we have: 



(2)°^ - ,(2)c 



(^ (0) (*)) E^^w) 2 ^^^)) 1 



(44) 



2 

^ 2)c w+3i?v? )c w = a, 2 4 2)c w-.9 2 (V 0) (m)) 4 2)c (*)-5 2 (V 0) w) E W )s (w? )s w +^ 13c wv>f )c (*) 



3=1 



J45} 



Note that the phases cancel out in the back-reaction 
on the inflaton field (in Eq. |4"4"]) . but not in any of the 
other equations. 

The physics which these equations describe is the fol- 
lowing: Since the system initially has no second order 
perturbations, it is the interaction of two first order 
modes whose wavenumbers add up to k that will source 
fluctuations of wavenumber k at second order. The sec- 
ond order perturbation at wavenumber k is affected by 
all first order modes. Hence, even though the effect of 
each individual first order mode is of the order e 2 , the 
large phase space of modes which contribute can lead to 
a large back-reaction effect [44{ . 

The above equations can luckily be reduced slightly. 
Since we have chosen not to excite the tf>^ s and 
modes, no second order sinusoidal fluctuations will arise, 
that is, <pffl s = ipr? s = for all n, and so there is no 
need to consider equations (|4*U| and (|%2"j) . Moreover, H 

can again be set to zero; and every term involving <pf 



or ipj in the interaction sum acting as a source in each 
equation can be set to zero. Also, as discussed above, 
the terms linear in the fields having as coefficient (fj) 2 
in equations (|4"0)l through (|4"5)l are negligible compared 
to the mass term of the fields, and thus can be dropped. 
However, recall that this conclusion was reached by im- 
posing a cutoff equal to the Hubble radius on the smallest 
scale excited to first order. Consequently, all sums in- 
volving interactions of first order modes acting as source 
terms for the second order perturbation modes can be 
performed up to n = ^ = L ^JY ~ 6 x 10 24 . 

Generalizing these equations from the 1+1- 
dimensional case to the higher-dimensional case 
this time is a bit more involved than simply replacing 
the wavenumber n by its vectorial expression m. In fact, 
the simple replacement works for every terms except for 
the interaction sum in each equation, which needs to be 
modified as follows: 



31,---,3d=0 



h + m 




[i{ 


ji = m 







r 3l---3 

jt, — ,3d=0 



hi,\ji- ni \ + h 1:n + ni , ji ^ ni 

2< 5fei,|ii-ni| +4 1 j 1 +ri 1 , h = "1 



ji,---,jd=0 



1 I Shun ,ji^0 

2 I 2^,0 , h = 



S k d ,\j d -n d \ + S k d ,j d +n d , 3d ^ n d 
2S k d ,\j d -n d \ + h d ,j d +n d , 3d = n d 



1 \h d ,j d , 3d + 

2 I 2S kd , , 3d = 



(46) 

(47) 
(48) 
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E ^wtf.xw 



26fci,o , ji = 



24 d ,o , 3d = 



(49) 



for the (0^ c , ipn", 4>^ C an< i tPtf' equations, respectively, in the case of d spatial dimensions. In particular, for d = 3, 

A(2)c 



,,(2) 



the sum for the (f>n equation (|4ip can be rewritten as: 



E W» 

i,i,fc=0 



+ 3 V ^ (1)c ^ (1)c + 3Vii (1)c ^ (1)c + V ib (1)c ib (1)c 



j,k=0 



fc=0 



+3 
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+3 I V ?i (1)c ^ (1)c + V iIj {1)c i!j (1)c 

+ ° \ Z^ ^ijk tp \i-n,\(j+n y )(k+n z ) + Z^ ^n»jfc^O(i+n !) )(fe+n« ! ) 



. i,j,k=0 



j,k=0 



Similarly, for the V'™ 2 ' )C equation (|43 
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j,fc=0 



(51) 



Finally, for the <j>^ c ipg equations (|44"|) and (|45|) . the sum at the end of each equation gets replaced by, respectively: 
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ojk j ' 2 
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fe=l 
LH 



,(l)c 
00A; ^OOfc 
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fe=l 



(52) 



From the above discussion, we see that the differential 
equations for the second order perturbation modes reduce 
to that of driven harmonic oscillators, with a constant 
period in the case of <jr- 2 '°, and with a time-dependent 
period in the case of •ip^ c . It is thus possible to solve 
them by making use of the Green's function method. 



For the 4>^ c equations, we use the causal Green's func- 



tion for a simple harmonic oscillator: 



G{s,t) 



—sin(m(t — s)) , for t > 
, for t < 



(53) 
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while for the tp^ c equation, we make use of the function: 
I £rM ath S(u ,q,m(t- s)) , fort>0 



G(s,t) 







fort<0 



(54) 

where u; and q are defined as above. 

In the following, we first analyze the back-reaction ef- 
fect on the background fields, and then move on to study 
the back-reaction of first order fluctuations on the per- 
turbation modes themselves [45| 



of radius which yields the result: 



V i 



Similarly, 



r'drdVL- 
i r 



(LX^v 2 ) 2 



48tt 



(57) 



B. Back-reaction on the cj> Background 



We first tackle the description of the mode 



whose 



evolution is dictated by equation (|44| generalised to 3 
spatial dimensions. This mode describes how second or- 
der perturbations of the (j) field will modify how the back- 
ground inflaton cjffl is oscillating around the minimum 
of its potential. From (|4"4")) and (|52p . we see that it is the 
growth of the amplitude of the first order perturbations 
in the ip field that will source this back-reaction. 

(2)c 

More precisely, solving the 4> Q equation using the 
Green's function method and making use of ([531 . we ob- 
tain the following solution: 



4 2)C (*) 



-1 



ds — sin(m(i — s))g 4> c cos(ms) 



Ojk 



1 1 (0 3 +jf:(*S 

LH 



fc=i 



(55) 



Note that the phases won't cause any major cancellation 
of the terms in the summation, since all terms are posi- 
tive. The expectation value of the sum can thus easilly be 
taken, which will simply add an additional factor of 1/2 
to the solution (that is, the expectation value of cos 2 9 
over one period). Thus, we can write: 

_„2, 



x cosh 



47T771 



(2L) J 



ds sin(m(i — s)) cos(ms) 



g4>c(l — cosms) 



- t 
8 ^ 



(i 2 + j 2 + fc 2 ) 1/2 



t x if + fc2 ) 1/2 



3 x ^ 1 

2^k 

k=l 



(56) 



We see that the time-dependent part of the sums factors 
out, which means that the summations can be performed 
independently of the time integral. This process of sum- 
ming is made much easier by replacing the sums with 
integrals, which is a resonnable approximation, since the 
range of the i, j, k variables is quite wide. Moreover, in- 
stead of integrating over a cube of length size ^f- in mo- 
mentum space, we integrate over one eighth of a sphere 
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i,fc= 



and 



i (f + k 2 y/ 2 
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k=l 



rdrdO- 



,1 



L\ x l 2 v 2 
4^3 



(58) 



dk— = In ( = — 

k V 2%/3vr 



(59) 



These sums are dominated by the small wavelength first 
order modes since these modes have the largest vol- 
ume in momentum space. Hence, since the momentum 
space is three-dimensional, it is clear that the first of 
the sums dominates over the two other. Indeed, for 
m = 10~ 6 (in Planck mass units), the first term con- 
tributes O(10 49 ), compared to O(10 25 ) and O(10) for the 
two others. The prefactor of the integral (|56l) therefore 

4 1 3 2 

becomes — 2H^3^m ' independently of L, and the integral 
itself can be performed numerically. 

As a stability and consistency check of our method, 
we study how the amplitude of the second order fluctu- 
ations in (f> affects the background field as we vary the 
value of the mass of the inflaton field in Figure [5] The 
first striking feature of these graphs is that, for all con- 
sidered values of m, the second order perturbation of the 
background mode becomes dominant before ip reaches 
the minimum of the potential, that is, within a short 
time relative to the time scale of the problem. 

Moreover, as m increases from 10 -7 to 10 -5 ' 5 , the frac- 
tion of the phase of 4>^ needed before the back-reaction 
term starts to dominate decreases. Equivalently, the 
value of mt when the amplitude of the second order per- 
turbations becomes of similar order as the background 
field (j)^ shrinks as the mass of the inflaton is increased. 

However, in none of the cases studied is the time 
elapsed before the back-reaction effects on the back- 
ground field becomes dominant shorter than 



(60) 



(which is 10 3 Planck times for the parameter values we 
have chosen), the other relevant time scale of the prob- 
lem under study which is the typical time scale of the 
motion in the tachyonic direction. Since the tachyonic 
field sources the instability at the end of inflation, its in- 
verse mass will set the time scale dominating the growth 
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FIG. 2: Comparison of the 4> background </>' ' (green) with the second order back-reaction ip (blue) as a function of Planck 
time, for values of the mass of the inflaton of m = 10 -5 ' 6 , 



10" b and m = 10"" 



is displayed in Planck units. Every plot ranges up to indentical phase of the background field (jr- ', that is, up to mt = n/2, 



in Planck units. cj> has units of mass, and 

(0) 



which allows for a comparison of the relative growth of 
becomes dominant over increases. 



compared to 



a(0) 



As m grows, the value of mt at which 



of second order perturbations. Varying to should thus 
affect the growth of second order perturbations in a very 
smooth and mild way, which is what Figure [5] shows. 



The 



mode is found to grow as 



an insta- 



bility that rapidly outruns the amplitude corresponding 
to the first order perturbation modes, which we found 
above to oscillate as ~ cos mi. A fit of the function 
—e af +b with b and c as fitting parameters, presented 
in Figure [3] for the case m = 10 -6 , is found to be in 
agreement with the analytical values over a time scale of 
?7i _1 . However, it is to be expected that during its very 
early evolution, the relative error between the fit and the 



numerical 



will be large. Indeed, 4>^ c (0) — 0, while 
the exponential fitting function can never reach zero. 

A plot of how the fitting parameters a and b vary with 
to is displayed in the bottom part of Figure [3J confirming 

(2)c 

the stability and smoothness of the evolution of the (j>Q 
mode under variation of the inflaton's mass. The param- 
eter a goes as ~ to, while the behavior of b was found 
to match a ~ 1/to function. Therefore, we find the fol- 
lowing fitting function for the growth of the second order 
back-reaction on the background: 



4 2)C 



-exp 



8.071e-5mr 



7.477e-6 
(to + 1.33e-S 



35.60 



(61) 



C. Back-reaction on the ip Background 



(2)c 

We now describe the mode Vm=0' wnose evolution is 
dictated by equation (|45|) generalised to 3 spatial dimen- 
sions. This mode describes how second order perturba- 
tions of the ip field will induce the growth of a spatially 
homogeneous contribution to the tachyonic field. At first 
order in perturbation theory the mode is obviously unex- 



cited, and so there is no non-trivial background to com- 
pare its growth against. Note that a second order ip 
background only develops since our choice of phases of 
the first order fluctuations breaks the ip — > —ip symmetry 
of the Lagrangian. Averaged over the entire universe, we 
would expect the symmetry to be restored. However, in 
a finite volume there is no reason why the phases should 
obey the symmetry. 

Combining the V'o equation and the Green's function 
(f54"l) . we obtain the following solution: 



ds M ath S(u} ,q,m(t - s))g (p c cos(tos) 

o m 

LH LH 



LH 
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12 5>, 



(lW(l)c 
OOfc ^00fc 
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(62) 



Substituting the solution for ^i 1 ^ and ipn C obtained in 
and ([3T|) . respectively, the sums to be performed 

(2)c 

turn out to be the same as above for (f> Q . However, 
there is one important difference: in the case of the back- 
reaction equation for the background <p field, the phases 
of the first order modes cancelled out. This will not be 
the case here. We will take this difference into account 
by performing the sums including random phases (see 
Appendix). Moreover, the M at h,S function appearing in 
the Green's function can be dealt with by means of the 
same kind of approximation performed for the M at hC 
function in ([2"T]) : 



M ath S - 



9 2 <ft -9 2 tt 



2m 2 ' 4m 2 



< sinh (1 — cos(z)) 

1 to 



(63) 
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■ a(m) 




8 07e-5m 









. b(m) 

-1it:?59-6/(m+1.3e-8>-35.6 




FIG. 3: Top panel: Fit of the amplitude of the second order 
back-reaction |<^o 2 ^ c | (dashed red line) with the function ~ 
e at +b (blue line), for m = 10 -6 . The fitting function is found 
to be in good agreement with the growth rate of I^q 2 ^! after 
mt — 0.3. The (f> background </>' ' (green) is also presented 
in order to provide the time scale of the growth of |<^q 2 ' c | and 
to provide a scale of its relative amplitude. Bottom left: 
Fitted value of a(m) as a function of m for values of m ranging 
from m = 1CT 6 ' 5 to m = 1CT 5 ' 6 (green points). The fit (in 
red) shows that the points obey a linear relationship over the 
studied range of m. Bottom right: Fitted value of b(m) as 
a function of m, over the same range. The fit (in red) exhibits 
that b(m) grows as — over the studied values of m. 



FIG. 4: Top panel: Fit of the amplitude of the second or- 
der back-reaction IV^'l (dashed red line) with the function 
~ e at2+b (blue line), for m = 10" 6 . The fitting function is 
found to be in good agreement with the growth rate of |</jq 2 ' c |. 
Bottom left: Fitted value of a(m) as a function of m for 
values of m ranging from m — 10™ 6 ' 5 to m = 10 -5 ' 6 . The 
fit (in red) shows that the points obey a linear relationship 
over the studied range of m. Bottom right: Fitted value of 
b(m) as a function of m, over the same range. The fit (in red) 
exhibits that b(m) grows as — over the studied values of m. 



is varied: 



Proceeding with these substitutions and performing 
the sums as a random walk, we get: 



33 87r 2 mL5 
x cos 2 (ms) cosh 



ds sinh 



3 (1 — cos(m(t — s))) 



cos ms 



m 



4 2)c 



-exp 4.139e-5mt 2 - '■ - 111.6 

(m + 4.781e-8) 

(65) 

This relation is displayed in Figure 2] for the case m = 
10~ 6 , where it is found to be in agreement with the ana- 
lytical values over a time scale of m , but as explained 
above, it deviates over smaller time scales. 



(64) 



D. Back-reaction on Long Wavelength <f> Modes 



which can now be integrated numerically. The result 
turns out to be quite similar to the one obtained for 
<Pq' : the second order background mode again grows 

.2 

as ~ — e , which is much faster than the first order 
perturbation modes, which we found above to grow as 
~ cosh(l — cos mt). However, due to the cancellation of 
terms in the sum induced by the random walk, their over- 
all amplitude is significantly smaller. In fact, even their 
rapid growth cannot compensate for this factor over the 
time scale set by the mass of </>. 

A fitting function very similar to the one found above 
for |^q 2 ' ) ' c | is found to match the growth of ,c | as m 



Wc now consider the case of main interest, i.e. the evo- 

(2)c 

lution of the large wavelength perturbation modes <f% 
(with ii small) at second order in perturbation theory. 
Their evolution is governed by equation (|41[) . with the 
interaction term replaced by (I50|) . However, once simpli- 
fied, the homogeneous version of these equations become 
exactly the homogeneous version of the equation we had 

(2)c 

for 0q , so that they can be solved by the same Green's 
function. Their solution thus reduces to the same inte- 
gral as before, and only the prefactor will differ. The 
main difference is that the cancellation of the phases ob- 

(2) 

served in the 4>o case is not present anymore, and sums 
over random phases need to be taken into account. The 
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expectation value over phases must thus be performed as 

(2) 

a random walk, as was done in the ip^ case. Terms with 
sums over 3 dimensions will again dominate over the ones 
with sums over only 1 or 2 dimensions. Thus, we are left 
with the task of evaluating: 
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(66) 



where each term can be summed independently as a dis- 
tinct sum (again with random phases). 

To evaluate these, we assume, for the moment, that all 
modes are in phase. Under this assumption, we compute 
the first and second sums separately, and note that the 
third and fourth will be bounded by the values of the 
two firsts. Since they turn out to differ by a negligible 
amount, we will be able to set them all to be equal to 
each other. We then restore the randomness of the phases 
to compute the sums under the approximation that they 
are equal their monopole term in a multipole expansion. 

Starting with the first sum, we note that it is bounded 
bellow by the sum (|57[) . which is in fact the monopole 
contribution to its value in a multipole expansion. Since 
we are considering long wavelength modes (m small) , set- 
ting the sum equal to its monopole contribution consti- 
tutes a first approximation of its value. However, we 
wish to know how good this approximation will be as a 
function of m. To answer this question, we compute the 
difference between the sum and the monopole term (|57[) , 
term by term: 



{i 2 + j 2 + k 2 ) 
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I7I-(* 2 +^ 2 ) 
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(67) 



This expression puts an upper bound on the term by 
term difference between the sum we want to evaluate and 
its monopole approximation (|5T[) . This difference can 
then be summed over the three indices i, j, k, in order 
to obtain an upper bound on the total difference. To do 
so, we again use integrals to approximate the sums, and 
perform them numerically. We obtain, for m = 1CP 6 in 




FIG. 5: Comparison on a log-log scale of the growth in ampli- 
tude of the <j> back-reaction on large-scale fluctuations (blue) 
with the growth of the first order ip perturbations, for m rang- 
ing from 1 to 10000 (red to yellow spectra), and for m = 10~ 6 . 
The evolution of the <f> background is also shown (topmost 
green curve) to provide the overall time scale of the evolution. 
Back-reaction becomes dominant within less than mt = 1 for 
all long wavelength modes considered (about 5 x 10 J Planck 
times). 
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and 



T 



didjdk 



\i 2 + j 2 + A: 2 ] 5/4 



-2.1 xlO 12 . (69) 



First, note that once the value of mis fixed, which mode 
exactly we choose to consider on the momentum shell of 
radius |m| will make at most a factor of \/3 difference in 
the result we just obtained. Hence we can take m = n x 
and n y = n z = for simplicity. 

Now, since the monopole term (|57| is O(10 49 ) for the 
parameter values which we are considering, we can be 
confident that the higher multipole corrections to the first 
sum are negligible as long as we choose O(ii) < 10 -12 . 
Since we are interested in long wavelength modes, we can 

take the first sum in (|66|) to be equal to — ■ 

If we now tackle the evaluation of the second sum in 
(1661) . we find that (I57|) is an upper bound on its value. 
To compute the difference between this upper bound and 
the actual sum, we proceed by following the exact same 
approximation scheme as above. In doing so, we find that 
the difference is bounded by: 
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FIG. 6: Comparison on a log-log scale of the growth in ampli- 
tude of the ip back-reaction on large-scale fluctuations (blue) 
with the growth of the first order ip perturbations, for m rang- 
ing from 1 to 10000 (red to yellow spectra), and for m = 10 -6 . 
The back-reaction never becomes dominant during the range 
of validity of our analytical approximations, that is, they do 
not become dominant before a time greater than mt = 3. 



which is once again negligible, provided that we choose 
0(m)<l(r 12 . 

Hence, we conclude that the four sums in (|66[) can 
be taken to be equal without making any significant er- 
ror. Now, performing the sums including random phases 
within this approximation, we find that the evolution of 
the back-reaction term for long wavelength modes (with 
O(m) < 10~ 12 ) is given by the solution: 



</4 2)c (t) = 
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'■(ftcvX 1 /* 



x cosh 



ds sin(m(i — s)) cos(ms) 
o 

cos ms) 



(71) 



Since this integral form of the solution is the same as the 
one we obtained for <p Q , the growth and m dependence 



of 



(for long wavelength modes) will be exactly as 



studied above for the back-reaction on the background. 
The only difference comes in the overall normalisation of 
this integral, which is now greatly reduced by the random 
phase, instead of being proportional to the volume of 
first order excited modes in phase space. Our conclusions 
concerning the consistency and stability as a function 
of m near the value m — 10~ 6 (in Planck units) thus 
still hold. In particular, we obtain the following fitting 
function: 



-exp 



8.071e-5mr 



7.477e-6 
(m + 1.33e- 



146.4 



Note that for long wavelength modes, it is independent 



of ii, and so the back-reaction term will have an identical 
growth for all modes with ii — 0(1O 12 ) . 

Comparing the growth of the second order contribu- 
tions to the amplitude of the long wavelength modes with 
the growth of the first order terms, we immediately see 

,2 

that, since the former grow as ~ — e , they will eventu- 
ally come to dominate over the corresponding first order 
terms. Indeed, the latter exhibit no instability since they 
simply oscillate as ~ cos mt. The precise time scale on 
which the second order terms become dominant can be 
found by comparing the the fitting function (fT2j) to ob- 
tain a semi- analytical estimate. 

Figure [5] shows the evolution of long wavelength fluc- 
tuations and the corresponding back-reaction term for 
various values of m, for m = 10 -6 . For this value of m, it 
takes between O(10 5 ) and O(10 6 ) Planck times for the 
amplitude of the second order long wavelength mode to 
become equal to the first order fluctuation of that mode. 
This time scale is comparable to the longest typical time 
scales of the problem, m~ , but is very long compared to 
the instability time scale in the tachyonic direction, m^. 

Since the only ii dependence of the intersection is 
through the normalization of the first order fluctuations 
modes in (|36[) (increasing by two orders of magnitude the 
momentum considered will reduce the amplitude of the 
first order perturbation by only one order of magnitude) , 
the time when the second order effects start to domi- 
nate over the first order perturbations will depend only 
slightly on the value of m considered, as long as only long 
wavelength modes are considered. 



E. Back-Reaction on Long Wavelength ip Modes 

Repeating the above analysis to compute the back- 
reaction on the ip perturbations, i.e. on ip^ c , for long 
wavelength modes we obtain (once again taking into ac- 
count that the phases are random): 
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which is simply 8ipQ ■ We thus obtain the following 
fitting function: 



V>i 2)c « -exp 4.139e-5mt 2 - , , "1 p , - 113.7 
v m y I (m + 4.781e-8) 

which is still independent of 11 for long wavelengths, i.e. 
for 11- O(10 12 ) . 

Again, the growth of second order modes will be much 
faster than the growth of the corresponding first order 
modes, since their instability only grows as cosh(l — 
cosm£). However, the amplitude of the back-reaction 
term is suppressed by the random phases, and so the 
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time interval At it takes before the second order term 
becomes more important than the first order term for 
long "wavelength modes is much larger than the time in- 
terval for the onset of the instability of the first order 
fluctuations. To find the exact crossing time, we need to 
find the intersection of (|75|) and (|57| . The result, shown 
in Figure [51 is O(10 6 ) Planck times for m — 10 -6 . This 
is much longer than the instability time scale for fluctu- 
ations of the waterfall field. Since once again the only 
m dependence of this crossing point is through the nor- 
malization of the first order modes, the chosen value of 
the long wavelength mode m will only mildly reduce this 
result as m is increased. 

Since At is much larger than the inverse mass of the 
waterfall field m^, we conclude that At is sufficiently 
large so that the back-reaction cannot shut off the growth 
of the entropy fluctuations before they have time to in- 
duce an appreciable contribution to the curvature per- 
turbations during the early stages of the tachyonic insta- 
bility. 



VI. CONCLUSIONS 

We have considered classical dynamics in a two field 
inflation model like hybrid inflation in which one field - <j> 
- is slowly rolling during inflation and a second field, the 
so-called "waterfall field" ip, develops a tachyonic insta- 
bility once <fi decreases below some critical value <f> c . We 
focus on the dynamics right after the instability develops, 
which corresponds to the initial stages of preheating at 
the end of inflation. 

We have provided a perturbative analytical analysis 
of the back-reaction effects of linear fluctuations on the 
background fields and on the perturbations themselves 
during the initial phase of the tachyonic instability which 
arises at the end of a period of hybrid inflation. We con- 
clude that the tachyonic growth of the short wavelength 
fluctuation modes leads - at second order in perturba- 
tion theory, the lowest order at which mode mixing ap- 
pears - to contributions which could in principle shut off 
the tachyonic growth of long wavelength entropy modes. 
However, in the case of random phases for the linear fluc- 
tuations considered here, we find that time scale when the 
back-reaction terms become important is much longer 
than the typical instability time scale of the tachyonic 
field. This process cannot shut off the instability of en- 
tropy modes before the latter have had time to become 
important. Note that there are other mechanisms which 
could dramatically weaken the strength of the instability 
of the long wavelength modes of the waterfall field, e.g. 
the decrease in amplitude of the linear fluctuations on 
super-Hubble scales during inflation. 

It is important to stress the crucial role which the as- 
sumption of random phases has played in our analysis. 
Had we worked with correlated phases, then the back- 
reaction effects would have been larger by a factor compa- 
rable to L 3 / 2 (in Planck units), and back- reaction would 



then be able to cut off the tachyonic instability of modes 
of the waterfall field long before they have time to de- 
velop. 

Our work complements the numerical studies of this 
model [U [22j ■ Our analytical approach has as a big ad- 
vantage that we can access the large dynamical range of 
scales which is required in order to draw conclusions for 
modes of cosmological interest today from dynamics tak- 
ing place on microphysical scales at the end of inflation, 
the scales which are studied in the numerical works. Dis- 
advantages of our approach based on a discrete Fourier 
space analysis are that it is purely perturbative and that 
we cannot study position space phenomena expected in 
the model such as the formation of topological defects. 

The reader may worry that the back-reaction effect 
which we find on the adiabatic mode of the fluctuations, 
the <f> fluctuations, might affect the predicted curvature 
fluctuations of the model. However, this is not the case. 
On scales larger than the Hubble radius, the curvature 
fluctuations are dominated by gravity. Microphysical ef- 
fects which occur at the time of reheating cannot affect 
the adiabatic modes of the curvature perturbations. This 
follows in full generality from the Traschen integral con- 
straints [39| which hold as long as we operate in the 
framework of General Relativity. In the case of scalar 
field models of inflation, it was shown in [ll[ that effects 
at preheating cannot affect the adiabatic mode of the 
curvature fluctuations. 

The reader should also keep in mind that our analysis 
does not concern other sources of entropy fluctuations in 
models like the one considered here. Specifically, if topo- 
logical defects are produced during preheating, then the 
active evolution of the defects between the time of forma- 
tion and the present time will lead to a large contribution 
to the entropy fluctuations. In fact, the specific model 
considered in this paper, where both <p and ip are single 
component real fields, would lead to domain walls and 
would in fact be ruled out because of the domain wall 
problem (40| . the over-abundance of energy in the walls. 
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Appendix: Sum over Random Phases 

In the evaluation of the sum of back-reaction terms, 
each term arises from an interaction between the two 
fields which each have a random phase and which, in 
general, do not cancel. We need to take the average 
over possible choices of phases not before studying the 
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effects on back-reaction, but only after. That is, the sum 
present in the expression for the background fields and 
their fluctuations must be performed for every possible 
choice of phases, and only then the expectation value of 
the obtained displacement distribution must be taken. 

Since we take the phases to be ramdom, this effect can 
be modeled by weighting randomly each amplitude, in 
the sum of the backreaction term, with 1 or -1. Calcu- 
lating the expectation value of the absolute value of this 
sum is similar to calculating the average displacement 
of a random walk in one dimension, so we will call this 
quantity a random walk summation. 

For a one dimensional random walk with N steps of 
size (5, the expectation value of the displacement (E\D\) 
is given by: 



hm m =s . 



(75) 



As the random walk summation of one shell, we thus 
obtain: 



47rr 2 1 1 



r 2tt(2 j L) 5 



1 



2 2tt(2£) 5 



(77) 



for each shell. This result is now independent of the 
radius, which allows us to take a random walk over the 
shells, that is a random walk of LH/tt steps with stepsize 
2-n(2L)' 1 ( n °te that just summing the above result over 
the shells is not the right way of proceeding, because the 
value calculated above could be + or — from one shell to 
the other, and doing so would thus greatly over-estimate 
the average result). We finally obtain: 



But in our case, this process is in fact slighly non- 
trivial, since the step size 5 is a function of the radius (r) 
in momentum space: S — > 5(r) — ~ ■ ^ or simplic- 

ity, we consider the random walk sum on the different 
shells of the sphere in momentun space independently, 
and we give each shell a thickeness of 1. The number of 
steps is then i the volume of the shell, that is: 



=i 



(76) 



LH 1 

7T 2tt(2L) 2 



87T 2 L 3 / 2 



(78) 



We also did a root mean square approximation to what 
the expectation value of the absolute value of the sum is, 
and we got the same expression within a numerical factor 
of order 1. 
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